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Abstract 

The algebraic structure of the ID Hubbard model is studied by means of the 
fermionic -R-operator approach. This approach treats the fermion models directly 
in the framework of the quantum inverse scattering method. Compared with the 
graded approach, this approach has several advantages. First, the global properties 
of the Hamiltonian are naturally reflected in the algebraic properties of the fermionic 
i?-operator. We want to note that this operator is a local operator acting on fermion 
Fock spaces. In particular, 50(4) symmetry and the invariance under the partial 
particle hole transformation are discussed. Second, we can construct a genuinely 
fermionic quantum transfer transfer matrix (QTM) in terms of the fermionic it- 
operator. Using the algebraic Bethe Ansatz for the Hubbard model, we diagonalize 
the fermionic QTM and discuss its properties. 
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1 Introduction 



Since the discovery of high T c superconductivity, strongly correlated electron models in low 
dimensionals have attracted much more interest. Among them the ID Hubbard model is 
one of the most interesting solvable models. The Hamiltonian consists of the hopping term 
and the on-site Coulomb interaction term, 

L L 1 1 

ft = - ( c ]+i* c j° + c ]aCj+i*) + U J2(n n - -)(n n - -). (1) 

j = l <T=U J' = l 

Lieb and Wu diagonalized the Hamiltonian [1] by using coordinate Bethe ansatz method 
under the periodic boundary condition (PBC). They found that the ground state at the 
finite correlation strength (U) shows property of the insulator in the half-filled case, and 
that of conductor in the other case. Some excited states are studied in [2, 3]. Takahashi 
introduced the string hypothesis to classify the states [4] and obtained the thermodynamic 
formulation The bulk properties, such as the specific heat and the magnetic susceptibility, 
were calculated using this formulation [5, 6]. 

In the frame work of the Inverse Scattering method [7, 8, 9, 10, 11], the Yang-Baxter 
equation plays an important role. Integrability can be discussed in the sense of an infinite 
family of conserved currents which are created from the transfer matrix. Then the transfer 
matrix is diagonalized by using algebraic Bethe ansatz. At the same time, this procedure 
can be applied to diagonalize the quantum transfer matrix which shows thermodynamic 
properties. Shastry introduced the Jordan- Wigner transformation [15] for the the Hubbard 
model to show the integrability of this model [12, 13, 14]. This corresponding spin model 
is called the coupled spin model (the Shastry model). Mapping to the spin model is a 
usual method since the ID quantum spin model can be easily reduced to the 2D classical 
spin model. Furthermore the quantum transfer matrix method was also applied to the 
coupled spin model [16]. The bulk properties, such as the specific heat and the magnetic 
susceptibility, were successfully calculated. They found that the spin-charge separation 
still remains even at the finite temperature. 

We have to note, however, that there exist several differences between the fermion 
system and the spin system. Because of the non-locality of this transformation, the PBC 
for the fermion system is not equivalent to the PBC for the corresponding spin system. 
Even when the system size is infinite, the correlation function, such as the one particle 
green function < cj^c^ > (a =f> !)> can n °t be expressed locally in the spin system and 
it is difficult to evaluate the thermodynamic properties by using quantum transfer matrix 
method. Thus another approach is required to investigate local properties. 

We introduced recently the fermionic /^-operator approach, where the genuinely fermionic| 
quantum transfer matrix can be constructed [21]. The fermionic -R-operator is constituted 
from fermion operators. The integrability is proved for the fermion models under the PBC 
[17, 18, 19] and the open boundary condition [20]. Here the super trace (Str) and the super 
transpose (st) for the fermion operators are introduced. For the spinless fermion model, the 
correlation length and the Uf oscillation are obtained at a finite temperature [21] through 
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the quantum transfer matrix method. This oscillation is a intrinsically fermionic property 
and the conformal field theory has already expected it. The correlation length also agrees 
with the finite temperature correction (T — > 0) [22]. For the Hubbard model, we showed 
the integrability using the the fermionic -R-operator in the previous study [19]. In this 
paper, we proceed this study and diagonalize the quantum transfer matrix algebraically to 
get the thermodynamic properties. The algebraic properties are also examined. 

We want to mention another attempt to treat fermion models directly, without mapping 
to the spin model. It is based on the graded Yang Baxter equation, where the quantum 
space is a fermion fock space and the auxiliary space is a graded vector space. The inte- 
grability [23, 24, 25] was first proved by Wadati et al. exactly for the Hubbard model. The 
SO (A) symmetry [26, 27, 28, 29, 30] of the transfer matrix is shown in [31, 32]. Ramos and 
Martins diagonalize the transfer matrix algebraically [33, 34], which agrees with the coor- 
dinate Bethe ansatz method [1]. Note that several differences exist between this approach 
and the fermionic -R-operator approach. First, the initial condition of the transfer matrix 
is not a shift operator in the graded approach and it is difficult to examine the energy 
momentum. Second, the fermionic -R-operator approach is more applicable when we apply 
the quantum transfer matrix method and diagonalize it algebraically. In this approach, 
both the quantum space and the auxiliary space are the fermion fock spaces. On the other 
hand, the quantum transfer matrix acts on the auxiliary space while the transfer matrix 
acts on the quantum space. Then both transfer matrices have a similar algebraic property 
and can be diagonalized almost in the same way. 

This paper is constituted as follows. In the next section, we review the fermionic 
.R-operator and the integrability of this model. In §3, the properties of the fermionic R- 
operator and the relation with global symmetry is clarified. In §4, we diagonalize the 
transfer matrix by means of algebraic BA approach. Then we generalize the eigenvalue 
of the vacuum state which is useful in diagonalizing the quantum transfer matrix in §5. 
The twisted boundary condition is also considered. In §5, we show that the monodromy 
matrix for the quantum transfer matrix is intertwined by the fermionic -R-operator just in 
the same as that for the transfer matrix. Then we diagonalize it algebraically and compare 
with [16]. The last section is devoted to concluding remarks. 

2 Fermionic i?-operator and Yang Baxter equation 

In this section we review the fermionic -R-operator approach for the ID Hubbard model 
(1). The integrability can be discussed by using the Yang-Baxter equation [19], 



&1 2 (U1, ^2)^13(^1, ^3)^23(^2, U 3 ) = % 3 {U2, U 3 )^i 13 (u 1 ,U 3 )^i 1 2(u 1 ,U2). 



(2.1) 



The monodromy operator 7 a (u) can be intertwined by the fermionic -R-operator, 



Rab(u a ,u b )7 a (u a )7 b (u b ) = 7 b (u b )3l a (u a )Jl ab (u a , u b ) , 
7 a (u) = R aL (u, 0)---XiK 0) 



(2.2) 
(2.3) 
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which leads to the commutativity of the transfer matrix t(u), 

[t(u),t(v)] = 0, r(u) = Sti a 7 a (u). (2.4) 

Here Str is defined as (B.8)(B.9). And the commuting family can be produced from the 
transfer matrix, 

[/ (m) > 7 (n) ]=0> = £_ {t{Q) ^) T ( U )}\ U=0 , (2.5) 

where 1^ is the Hamiltonian. The fermionic /^-operator which satisfies (2.1) [35, 36, 37, 19] 
is given by 

3% 2 (u,v) = R ( £(u - v)$$(u - v) + ~ V J tanh(%) - h(v)) 

x ^ 2 \u + v)^ 2 \u + v)(2n^ - l)(2n u - 1), (2.6) 

$$(u) =a(u)(-n lr7 n 2r7 + (1 - n la )(l - n 2tr )) - b(u)(n la (l - n 2a ) + (1 - n la )n 2a ) 

+ c{u){c{ a c 2a + cl a c la ), (<T=tl) (2.7) 

sinh2/i(-u) [/ 

where a(u) = cosw, o(tt) = sinu, c(u) = 1, = — , (2.8) 

which is slightly different from the one in [19] because of the different JW trans. Note 
that the initial condition for this /^-operator (u=v=0) gives the permutation operator (y), 
which means that the initial condition for the transfer matrix r(0) is a left shift operator, 

^jkCja c/c,j!P ' jki '■PjkCkij CjcT^jk a nd Tj^ 1, (2-9) 
r(0)c ja = c j+la r(0), with PBC (c L+j = Cj ), (2.10) 

The Hamiltonian (1) is expressed as 

d L d 

H = r{Q)- l -r{u)\ u= , = S^+i^j+i^OlU (2.11) 

j=i 

This is the fermionic i?-operator approach for the integrability. We see that (2.2) (2.3) 
can be generalized [19] to 7 a (u) = Jl a L(u,u ) ■ ■ - 3l a i(u, u ) where Rjk(u , u ) = O 5 ^. Then 
the generalized Hamiltonian is obtained. It was generally believed that we demand the 
deference property Jl(u, v) = "Jl(u — v) to get the boost operator. The Hubbard model does 
not have this property. Recently J. Links et al. found the boost operator by using this 
generalized formula [38]. 



3 50(4) Symmetry and Partial Particle-Hole Trans- 
formation 

In this section we discuss the two important symmetries of the Hubbard model by means 
of the fermionic i?-operator. The first one is the 5*0(4) symmetry which consists of the 
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spin- SU (2) and the charge- SU (2) (^-pairing SU(2)). The generators of the spin-5'C/(2) 
and the charge-S'L r (2) are given by 

L L , L 



rf 



j=i j=i j=i 

(3.1) 

1 L 1 L 



2 "Jim " ' 2i 

i=i j=i 



1 L 

^ = ^E^t + ^r 1 )' (3-2) 

These generators constitute the Lie algebra of 5*0(4), i.e., so(4) = sw(2) © sw(2). It is well 
known that the six generators (3.1) and (3.2) commute with the Hubbard Hamiltonian (1), 

[H,S a ] = [H, V a ]=0, (a = x,y,z) (3.3) 

if we assume the even number of sites and the PBC. 

The second one is the so called partial particle hole transformation, 

c n ^(-iyc] L , U^-U. (3.4) 

It is also well known that the Hubbard Hamiltonian is invariant under the partial particle- 
hole transformation, while the generators S a and i] a are exchanged to each other. 

We discuss the above two symmetries in terms of the fermionic /^-operator. Here 
we omit some calculations, which we will show in Appendix A. The Lie algebra so(4) 
symmetry of the fermionic i?-operator was already discussed in [19] , which is represented 
by the relations, 

[Jl jk (u,v),S? + S%] =0, (a = x,y,z) (3.5) 
[% jk (u,v),r j z j +ri*\ =0, (3.6) 
{% jk (u,v),ri?-r$} = 0. (a = x,y) (3.7) 

Here we have introduced the local generators as 

1 1 1 

rf! = \( c )Ai + c n c n)> Vj = jX c ]Ai + c n c n)> = ^Kt + n jl - 1). 



The relations (3.5)-(3.7) can be extended to the those for the transfer matrix, from which 
it is possible to prove the so(4) symmetry of the local conserved currents. 
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We shall generalize the symmetry relations (3.5)-(3.7) into exponentiated forms, which 
corresponds to the generalization of the Lie algebra so(4) symmetry into Lie group 5*0(4) 
symmetry. By direct calculation, we have found the following identities, 

[R jk (u, v),exp(i9Sf) exp(i0S£)] =0, (a = x, y, z) (3.8) 
pR jfc («, v), exp(i^|) exp(i^)] = 0, (3.9) 
Rjk(u, v) exp(i6r]J) exp(-i^) = exp(-i^) exp(i6>0% fe (w, v), (a = x, y) 

(3.10) 

Here 9 is an arbitrary (real) parameter. From the symmetry relations (3.8)-(3.10) for the 
fermionic -R-operator, we can derive the following identity for the transfer matrix, (see 
Appendix A) 

exp(— i9S a )r(u, u ) exp(i9S a ) = t(u, u ), (a — x,y,z) (3-H) 
exp(-i9r) z )T(u, u ) exp(i#rf ) = t(u.u ), (3.12) 
exp(i^ Q ) r(-u, u ) exp(i9r) a ) = r(u, u ). (a — x,y) (3.13) 

When we think about the local conserved currents (J*™)) as was discussed in (2.5), we 
consider the product r _1 (-u , uq)t(u, uq) and it satisfies the relation 



c 



{9x t-\uq, u )t{u, u )e iex = r _1 (uo, u )t(u, u ), (3.14) 



where X is any generator of 50(4), i.e., S a or i] a (a = x, y, z). Namely the local conserved 
currents which contain the Hamiltonian are therefore invariant under the 50(4) rotation. 

Next we consider the partial particle- hole transformation (3.4). The relevant relations 
for the fermionic -R-operator are found to be 

% jk (u, v- U)VjV k = VkVjXjkiu, v- -U), (3.15) 
% jk (u, v- U)VjV k = V k Vj3l jk (u, v- -U), (3.16) 
where Vj = (1 - 2n jT )(c} i + Cji ), Vj = i(l - 2n jT )(c J t i - c u ) = i(2n n - l)V r (3.17) 

Here we write the [/-dependence of the -R-operator explicitly. Note that Vj and Vj are 
Grassmann odd operators and Vj = Vj = 1. The operators Vj and Vj induce the following 
transformation for the fermion operators, 

VjCj^Vj = c jh VjC^Vj = c) v \'j<- i \'j = -c) v VjC^Vj = -c^. (3.18) 
Now let us define the operator V by 

L/2 

V = l[ V 2k V 2k - 1 = V L V L . X ■ ■ ■ V 2 V U (3.19) 
k=i 
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From the property (3.18), one can see that the transformation for an arbitrary operator 
X({c ja },U) 

X({c 3a }, U) — V~ l X({c 3u }, -U)V (3.20) 

is nothing but the partial particle-hole transformation (3.4). From the properties of the 
fermionic -R-operator (3.15)(3.16) the product r _1 (-u , u )r(u, u ) is invariant under the 
partial particle-hole transformation, (see Appendix A) 

V-W-^Uo, u )t(u, u ; -U)V = r-\u , u )t(u, u ; U). (3.21) 

Then we can conclude that all the local conserved currents are invariant under the partial 
particle- hole transformation (3.4). 

4 Algebraic Bethe Ansatz 

In this section, we propose the brief note of the algebraic bethe ansatz. As Ramos and 
Martins did in the graded case [33], we construct the eigenstates algebraically although a 
little difference exists. The eigenvalue of the transfer matrix and the nested bethe ansatz 
equation are obtained. Then we apply the generalized eigenvalue of the vacuum state, 
which we will use in §5. Twisted boundary condition is also discussed. 
We will diagonalize the transfer matrix (B.8)(B.9), 

t{u) = Str a 7 a (u) = D 11 (u) + D 22 (u) - A u (u) - A 22 (u). (4.1) 

The eigenstate is constructed algebraically. We have five creation operators (C 2 i, C 22 , D 21 , 
Bu and B 2 i) defined in (B.7). By using three creation operators (C 2 i, C 22 , D 2 \) and taking 
liner combination of these products, we define the n-particle eigenstate by 

|$n(Ai, • • • , A n )> = 4(Ai, • • • , A n ).#|0>, (4.2) 
<S„(Ai, • • • , A n ) = C{X X ) ® $ n -i(A 2 , • • • , A n ) 

+ £? =2 [|*<g> ^2i(Ai)$„_ 2 (A 2 , ■■ ■ - ■■ , A„) J D 22 (A,)].^ ) 1 (A 1 , • • • , A n ) (4.3) 

— * — * 

where C = (C 2 i,C 22 ) and £ = (0,1,-1,0). C 2 \ and C 22 are 1-particle creation operators 
and D 2 i is 2-particle creation operator. §j is a operator defined recursively, 

<?j-i(Ai, • • • , A n ) = j lj j g^ 2 (\!, • • • , Xj, Xj-i, • • • , A„).fj_ij(Aj_i, Xj), (4.4) 
(n) n X] _ /(Ai, A 2 ) a + (A fc , A 2 ) 

The vacuum state (|0)) is a fermion Fock space defined as, 

L 

|0> = ni°% where c jCT |00>j = (a =U), (4-6) 

3=1 
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The main point of this algebraic construction is its exchange property , 

$ n (Ai, • • • , A n ) = a 4(Ai, • • • , \j, Xj-i, ■■■ , XJ.rj-uiXj-i, A,), (4.7) 

a + (A j _i, Xj) 

which simplifies this construction. After a long patient calculations (see Appendix B) using 
the algebraic relations (B.10)-(B.19), we get the eigenvalue of transfer matrix 

Hu, {M ) = (o+ („))^ l( -^M) + ( c+ („))^ l( -£j^>) 

-( 6 -w) I ( n -(-^ )n£ '^ +n "- ( -^ )nE ^ } ' (4 ' 8) 

where a ± (?x) = a ± (u,0), b ± (u) = b ± (u,0) and ^(u) = (^(u, 0). And the Bethe Ansatz 
equations are, 

To see the agreement with the coordinate one [1], we examine the energy spectrum and 
the momentum, after the reparametrization e 2h( ~ x ^ / tanAj = e ikj , 

E= ^logA(«,{A,})| tl=0 = -2E^ =1 cos% + ^(L-2n) (4.10) 

P = ilog(A(u = 0, {A,})) = -EjUkj. (4.11) 

Now we compare with the Ramos and Martins' results. As we mentioned in §2, the initial 
condition of the transfer matrix doesn't give the shift operator in the graded case. Then the 
energy momentum is not given as (4.11). Beside the agreement of (4.9)(4.10), (4.8) differs 
from their results by the overall factor. As is discussed in [33], the completeness can be 
discussed in the similar way [39, 40]. Here the regular bethe ansats are the highest- weight 
states of the 5*0(4) symmetry. The other states are obtained by the lowering operators. 
Next we generalize the eigenvalue of the vacuum state in the following way, 

D u (u)\0 >= du(u)\0 >, D 22 (u)\0>= d 22 (u)\0>, (4.12) 
An(u)|0 >= a n a xx (u)\0 >, A 22 (u)\0 >= a 22 a xx (u)\0 >, (4.13) 

where an and a 22 are constant. The eigenvalue of transfer matrix A(u, {A/}) and the Bethe 
Ansatz equations are, 

Mu, {A,}) = MU -p^) + dn{u) n U - b ^M) 

(4.14) 

(^) l = an n Eir -i-, n; U %„A,) = ^n & ,§^. (4.45) 

\axx{Aj)' K^uXj) a 22 ^ b(iJ, k , Hi) 
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Then we think about the twisted boundary condition, 

cl+it = e Wl c lh c L+ll = e W2 c n , (4.16) 

where the transfer matrix is defined as 

T twist (u) = Str a 0^ wist (w), (4.17) 
0^ wist ( M ) = V7 a {u) where V = e^^^ 11 ^ 1 ^^-^). ( 4 . 18 ) 

Because of the properties of the fermionic i?-operator (3.8), the global Yang Baxter relation 
(2.2) is satisfied with this monodromy matrix < J twlst although the eigenvalues of its elements 
for the vacuum state change to 

d 22 ( u ) = a + {u) L , d u (u) = e i ^c + (M) L , (4.19) 
a n = e- ie \ a 22 = e~ w \ a xx {u) = b~ (u) L . (4.20) 

Here 9i = 9 2 = gives the PBC (cj +Lcr = Cj a ) and 6i = 6 2 = -n gives the anti periodic 
boundary condition (cj + L a = —Cj C ). 

5 Algebraic Bethe Ansatz for the Quantum Transfer 
Matrix 

In this section, we will diagonalize the quantum transfer matrix in the same way as §4. 
First we compare the quantum transfer matrix with the usual transfer matrix. Then we 
make use of the similar algebraic structure to diagonalize it. 
In the row to row case where the transfer matrix, 

t(u) =Str o 0l o£ («,0)--.a o i(u,0) =r(0)(l+uH + O(u 2 )), (5.1) 

is diagonalized, we use the Yang Baxter equation to intertwine with respect to auxiliary 
space, where 1 — * a, 2 — > b and 3 — > j in (2.1), 

Rab(Ua,Ub)Raj{Ua,Uj)%j(Ub,Uj) = %j(u b , Uj)3l aj (u a , Uj)Jl ab (u a , U b ) . (5.2) 

But in the case of the quantum transfer matrix, it operates on the auxiliary space and 
the monodromy matrix should be intertwined with respect to the quantum space. We use 
the two equations, 

$>jk(Uj, U k )$. ak (u a , U k )Jl aj (u a , Uj) = %aj(u a , Uj)3{ ak (u a , U k )Jl jk (u j: U k ), (5.3) 
%jk(Uj, U k )Jl ak (u k , Ua^ajiUj, U a ) = R a j{Uj, U a )5i ak (u k , U a )Jt jk (Uj, U k ), (5.4) 
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which are equivalent to (2.1). In (5.3), l.h.s. and r.h.s are reversed and 1 — > a, 2 — > j and 
3 — > k. In (5.4), we take the super transpose (st) with respect to 3 and 1 — ^ j, 2 — > A; and 
3 — > a, where the super transpose (st) is defined as st a = st\ st£ and 

(Xi/w + X 2 cL + X 3 c a(T + X 4 (l - n a(T )) st « 

= Xin a(T -X 2 c a(J + X 3 cl a + X 4 (1 -nj, (a =f, |) (5.5) 

where a operators in space a is expressed explicitly. ^k a j(uj,u a ) can be expressed in the 
following way, 

Raj(Uj,U a ) = $.%(Uj,U a ) (5.6) 

= Sfltuj - u a )^{ Uj - u a ) + C0S( , Mj 7 Ua \ tanh(fr(%) - %„)) 



COs(Wj + M a ) 



x ft^ + Ma )^a)( % + Ma )( 2%T - l)(2r^ - 1), (5.7) 



ftgV) = Wjf^u) = a(u)(-n Ja n aa + (1 - n ja )(l - n aa )) 



- b(u)(n ja (l - n aa ) + (1 - n jc )n aa ) + c(u)(c] a cl a + c ja c aa ). (5.8) 
Because of the properties for ft^(O), 

-^■aj+i( u )\u=oKj = HJ j+1 ?l j+1 9l j: (5.9) 

^'aj+l^aj C j+l'y = C j^aj+l^aj-> (5.10) 

where ^ = ^(0), HjJ +1 = -(c+ la c ia + c£c j+lCT ), (<r =f, I) (5.11) 

we have the properties, 

?(u) = Str a ft ai (0, «) • • • ft al (0, u) = (1 - uH + 0(« 2 ))f (0), (5.12) 
f(0)c j+la = c ja f(0). (5.13) 

Here f (0) is a right shift operator. Then 

r (0)r(0) = 1, (5.14) 
t(-u)t(u) = l-2uH + 0(u 2 ). (5.15) 

Next we evaluate the partition function by using ft and ft, 

Tre-^^ lim Tr(r(- M )f (w))^ 2 (5.16) 

A 1 "— >oo 

N/2 

= ^ OQ Tr II Str « 2m ,a 2m _ 1 [fta 2m L(-M,0)---ft a2m i(-M, 0) 
m=l 

xi^KOj-V.iKO)] (5.17) 

= Jim StrnTr,0f™K0), (« = |) (5.18) 
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where N is called Trotter number which we to be even and 

7f™(u,v) = 3l aN j( it, v)3l aN1 j (t>, u) ■ ■ •% a j(-u,v)% ll j(v,u), (5.19) 
r Q™(u,v) = Tr 3 7f™(u,v). (5.20) 

Here r QTM (u,v) is the QTM. From (5.3)(5.4), this monodromy operator is intertwined as 

%jk(uj, u k )7%™(u, u k )7f™(u, Uj) = 7f™(u, Uj)T^™{u, u k )Ji jk (u j: u k ). (5.21) 

After gauge transformation (B.3), we have the property for the fermionic -R-operator 
$lj k (uj,u k ) = ^kji.— u ki — u j)- Then we have 

%, jk (uj,u k )7j™(u, -Uj)7^™(u, -u k ) = 7%™(u, -u k )lf™(u, -u^j^u^Uk). 

(5.22) 

We see the similar relation for the usual transfer matrix (2.2) (2.4). And the QTM is 
diagonalized in the same way except the generalization of the vacuum eigenvalue. Next we 
introduce the magnetic field (H = 2h) and the chemical potential (//), 

L 

H cxt = - X^Kt - n 3i) - M^t + n jit )}, (5.23) 

and evaluate the QTM algebraically. Tr^ in (5.20) is nothing but the anti periodic boundary 
condition (9i = 9 2 = ir in (4.16) and V apbc in (4.17) ) and the introduction of the magnetic 
field and the chemical potential is also applied to (4.17) as V m ' c , 

V = v apbc V m ' c , (5.24) 

^apbc _ e -vri(n aT +n ai ) ym,c _ e -/3jt(n aT +n al )+/3h(n aT -n al ) ^ 

And the vacuum state (|fi) aux ) in the auxiliary Fock space is defined as 

N/2 

l^) aux = II 1°°)^ I 11 ) — ' where |H) a = cl r cij00> a (5.26) 

m=l 

and |00) a is given in (4.6). Then C 2 i, C 22 and D 21 remain to be creation operators. Then 
the generalized eigenvalue for the vacuum state is given by 

d 22 {v) = (a + (-u,-v)c + (-v,u)) N / 2 , d n (v) = e- 2 ^(c + (-u,-v)a + (-v,u)) N / 2 , 

(5.27) 

a u = e ^ +h \ a 22 = e ^~ h \ a xx (v) = (-b + (-u, -v)b~(-v, u)) N ' 2 . (5.28) 
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And the eigenvalue and the bethe ansatz equations (4.14) (4.15) are 




b{n h u) 



1 



(5.31) 



(5.29) 



(5.30) 



When the system size is infinite [L — > oo), we see that the largest eigenvalue of the QTM 
with (y = 0) gives the partition function because (5.18) can be expressed as a summation 
over the eigenstates (Str), 



where pi is a parity of n. The largest eigenvalue (Ai) is given by the sector (n — N 
and m = N/2) and this agrees with [16] after the partial particle hole transformation. 
Furthermore we can deal with the intrinsicly fermionic properties which can not be expected 
in [16]. The asymptotics of the correlation functions are given by another eigenvalue. For 
example, the eigenvalue of (n — N — 1 and m = N/2 — 1) sector gives the one particle 
Green function < Cj a Cj a >. (see [21]) 



We studied the algebraic structure for the ID Hubbard model in terms of the fermionic 
i?-operator. This approach is a powerful method when we treat fermion system directly. 
The quantum transfer matrix (QTM) is naturally constructed. Then it can be diagonalized 
just in the same way as that for the usual transfer matrix. It has been a difficult problem 
and mapping to the spin system was a usual method, where intrinsicly fermion properties 
can not be observed. Our result can be applied to evaluate any correlation function in 
addition to the bulk properties. 

In this paper, we review the fermionic i?-operator approach for the ID Hubbard model 
and show the supertrace and supertranspose for the fermion operator. Then the algebraic 
aspect of this approach are discussed. The properties of the i?-operator itself reflect the 
global properties of the model. Then we construct the eigenstate algebraically and diag- 
onalize the transfer matrix. There appears naturally the r-structure which leads to the 



Tre"^ 11 = A[ + (-1) P2 A^ + (-1) P3 A^ + • • • , 



(5.32) 



6 Concluding Remarks 
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nested bethe ansatz equations. The QTM is made from the fermionic i?-operator and the 
supertrace of it. By using the generalized algebraic bethe ansatz equations, we diagonalize 
the quantum transfer matrix (QTM). 

For the future work, we will investigate more details about the solutions of this Bethe 
ansatz equations for the QTM to obtain the thermodynamic properties of this model. 
Here the auxiliary space is finite and the Trotter limit (N — > oo) is desirable to obtain the 
exact result. For the bulk properties, Juttner et al. [16] obtained the nonlinear integral 
equations (NLIEs) to evaluate the eigenvalue of the QTM, where the infinite space can 
be easily considered. We expect that the NLIEs for the another eigenvalue can be the 
small deformation of the original one [21]. The one particle Green function c^Cka and the 
spin correlation function S^~S^ are our first problems, where we expect no singularity with 
respect to the temperature. The spin correlation function SjS%. is our next problem. We 
expect it to be more complicated. When the Coulomb interaction is strong in the half-filed 
case, it behave as the XXX Heisenberg chain, where the crossover phenomena in SjS% is 
observed [41]. 



The authors are grateful to M. Shiroishi, A. Kluemper and M. Wadati for useful discussions 
and careful reading of the manuscript. This work is supported by Grant-in-Aid for JSPJ 
Fellows from the Ministry of Education, Science, Sports and Culture of Japan. 



Here we will show some calculations which are omitted in §3. The global properties are 



obtained from the local properties of the fermionic i?-operator. As for the 5*0(4) symmetry, 
it is from (3.8)-(3.10) to (3.11)-(3.13). (3.11) is proved from (3.8), 
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Appendix A; SO (4) symmetry 



L 



exp(-i9S a )T(u,u ) exp(i9S a ) = Str, 



a 



L 



Str, 



a 



Y[e ws «R aj (u,u )e 



t(u,Uq). 



(A.l) 
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In the same way, (3.12) is obtained from (3.9). And (3.13) is from (3.10), 
exp(i9r] a )T(u, u ) exp(i6r) a ) 

L/2 

= Str a H e ie ^X a:2k (u, uoje^e-^Vi^- u )e~ id ^ 

k=l 

L/2 

= Str a H e ier, ^ a , 2k (u, ^e^e-^ViK u )e~ i9 ^ 
k=i 

= r(u,u ). (A.2) 

When we consider r _1 (-u , -u )r(-u, u ), it has 5*0(4) symmetry (3.14). 

As for the invariance under the partial particle hole transformation, we introduce V 
defined by 

L/2 

V = H V 2k V 2k ^ = VlVl-! ■ ■ ■ V 2 V 1 . (A.3) 
fc=i 

Now using (3.15)(3.16), we obtain 

L/2 

V" -1 r(«,u ; -U)V = Str a JJ V^^fcK^o; -U)V 2k V 2k _ 1 Jl i a:2k _ 1 (u,uo)V 2k _ 1 

k=i 

L/2 

= Str a i[(-VX,2 k (u,u 0] U)V a )(-VX,2 k -i(u:U )V a ) 
k=l 

= -r(u,u ;U). (A.4) 

Here in the last equality we have used a property of the supertrace 

Str a {V a XV a } = -Str a {X} . (A.5) 

for an arbitrary operator X. The transformation law (A.4) means that the transfer matrix 
itself is not invariant under the partial particle-hole transformation (3.4). When we consider 
the product t~ 1 (u ,uq)t(u,uo), as before, it is invariant (3.21) under this transformation. 

Appendix B; Fermionic i?-operator and Algebraic 
relations among the elements of the Monodromy 
matrix 

Here we show the simplified fermionic i?-operator and the algebraic relations among the 
elements of the monodromy matrix 7 which are necessary when diagonalizing the transfer 
matrix. The relations are obtained from the global Yang Baxter relation, 

Rab(u a ,u b )7 a (u a )7 b (u b ) = 7 b (u b )7 a (u a )3l ab (u a , u b ) (B.l) 
7 a (u a ) = 5l aN (u a , 0) • • • ftaiK, 0). (B.2) 
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Then some calculations that are omitted in §4 are supplemented. 
We gauge transformed the fermionic i?-operator (2.6), 

Ia(u a )Ib{Ub)Rab(Ua, U b )I a (u a )~ l 7 6 (w&) _1 -> %b(u a , U b ), (B.3) 

where 

I a (u a ) = e M«a)(2n oT -l)(2n 6i -l) ) ( R4 ) 

to the simple form, 

fti 2 (wi, w 2 ) = a + (Mi,M 2 )(ni T n u n 2T n2i + (1 - ni T )(l - n u )(l - n 2T )(l - n 2i )) 

- a~(ui, « 2 )(ni T (l - n u )n 2T (l - n 2i ) + (1 - ni T )n u (l - n 2T )n 2i ) 
+ c + (M 1 ,M 2 )(ni T n u (l -n 2T )(l - n 2i ) + (1 - n 1T )(l - n u )n 2T n 2i ) 
+ c _ (ui,u 2 )((l -ni T )n u n 2T (l - n 2i ) +ni T (l -ra u )(l - n 2T )n 2i ) 
+ b + (u l7 u 2 )((l - m T )(l - n u ) - ni T n u ) 

x (n 2T (l - n 2i ) + (1 - n 2T )n 2i ) 
+ 6~(ui, tt2)(ni T (l - n u ) + (1 - n n )n n ) 
x ((1 - n 2T )(l - n n ) - n 21 n 2l ) 

- d + '(u 1 ,u 2 )(c\^c\ l c 2 ^C2i + Ci T Cij.4 T 4|) 
+ d~{u u u 2 ){c\ ] c n c 2 ^c\ l + Ci T C U 4 T c 2|) 

+ e(Mi,M 2 )((c{ T c 2T + 4 T Ci T )((l -n H )(l - n n ) - n n n 2l ) 

+ (c f u c 2i + 4| c u)(( 1 - ^itX 1 - ^2t) - rai T n 2T )) 
+ /(mi,M2)((4 t c 2T + 4 T Ci T )(n u (l - n 2l ) + (1 - n u )n 2i ) 

+ ( c i| c 2| + 4| c u)(^it( 1 ~ n 2t) + (1 - ni T )n 2T )), (B.5) 



where 



a ± (uj,Uk) = cos 2 (uj — Ufc) 



1 ± tanh {h(uj) — h(uk)} 



cos 



COS Mi 



Uj + u k ) 
u k ) 



b (uj, Uk) = — sm(uj — Uk) cos(uj — Uk) 
= — sin(-Uj — Uk) cos(uj — Uk) 



1 ± tanh {h(uj) — h(uk)} 
1 ± tanh {h(uj) + h(uk)} 



sin(uj + u k ) 
sin(uj — Uk) 
u k ) 



COS Mj 



COS Mi 



u k). 



(^(uj, u k ) = sin 2 (uj - u k ) 



1 ± tanh {h(uj) + h(uk)} 



sin(uj + Uk) 
sin(uj — Uk) 



d ± (uj, Uk) = 1 ± tanh {h(uj) — h(uk)} 



cosjuj - u k ) 

COs(Uj + Uk)' 



= 1 ± tanh \h(uj) + h(u k )} — H 

sm{Uj + u k ) 



e(uj,u k ) 



cos ru.; 



u k ) 



cosh {h(uj) — h{uk)Y 



f(Uj,U k ) = 



sin(uj — Uk) 



cosh{h(uj) + h(uk)} 



(B.6) 
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We define the Monodromy operator as follows, 

7(u) = Dun^rii + D 22 (l - ^tX 1 ~~ n i) + ^n n t(l _ n i) + ^22(1 ~ n \) n i 
+ Dnc\c\ + D 2 \c^ci - A 12 c\c l + A 2 iC T c[ 
+ Cnn T cj + Ci 2 c|nj - C 2 ic-f(l - raj + 622(1 - ra T )q 

+ Bnn^ci - B 12 c\(l - ni) + B 2 \c^ni + B 22 (l - ra T )cj. (B.7) 

Note that the sign is a little tricky as to get the simple r-structure which is the most 
interesting point of this story. As the definition of the supertrace Str, the transfer matrix 
is expressed as, 

r(u) = Sti a 7 a (u) (B.8) 
= D u (u) + D 22 (u) - A n (u) - A 22 (u) (B.9) 

Relations among An ■ ■ ■ D 22 are 

C 2a {u)C 2b {u>) = a ^ U 'h f d (u,u')C 2c (u')C 2d (u) 
a + (u, u ) 

+ f ^ U '\ u{D 2l {u)D 22 {u>) - D 21 (u')D 22 (u)), (B.10) 
c + (u, u') 

A ab {u)C 2c {u') = - U '\ r b e c d (u, u')C 2e (u')A ad (u) + ^%C 2b {u)A ac {u') 
(u,u') (u,u') 

f(u,u') ,„ 1 » ... e(u, u') ^ . . „ ... 
+ > {B al (u)D 22 (u') - v ' > D 21 {u)B a2 {u') 

c + (u,u') b~(u,u') 

+ ^^D 21 (u')B a2 (u)}, (B.ll) 
(u, u') 

D 11 (u)C 2a (u') = _^l^C7 2a («0Ai(«) - ^^AnK)C la H 

+ ^d DAu)CUu) 

+ ■ / . ( "-"' ) ) {B 11 ( M )^ a («') - S 21 («)^ 1o («')}, (B.12) 
^(uOCaaM = f^ C^MJW) " pT-^ C 2a(u r )D 22 (u), (B.13) 
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A ab (u)D 21 (u') = -D 2l (u')A ab (u) 
e(u, u' 



, UC 2b (u')B al (u) - C 2b (u)B al (u')), 
D LX {»)D 2X („') = - d LS^ D2l{u)Dll{u >) + ^p^l D2l{u ' )Dll{u) 



+ 



c + (u, u') 

J^lf ab B al (u)B bl (u'), 
c + (u, u') 



c + (u, u') 



D 22 (u')D 21 (u) = - d -S^D 21 {u')D 22 {u) + a L^^D 21 {u)D 22 {u') 
■ "" ' c + {u, u') 



f(u,u' 



c + (u, u') 



y:t,abC 2 a(u')C 2b (u), 



B nl (n)C 2 M = ^^C Vt (u')B dl (u) 



+ 



b (u,u') 
e(u, u') 



b~(u, u 



J -(D 21 (u')A ab (u) - D 21 (u)A ab (u')), 



B a2 {u)C 2b {u>) = J^p^C 2b (u')B a2 (u) 

[UjU J 



e(u, u') 



b (u,u') 



J -{D 22 (u')A ab (u) - D 22 (u)A ab (u')), 



C la (u)C 2b (u>) = - d -^lc 2a {u)C lb {u') 



+ lp^MMA 2h {u>) - A 2a {u)A lb {u')) 
+ ?-^x%(u,u')C 2c (u>)C ld (u) 



/(«,«') 



c + (u, u') 



^U{D 22 {u')D n (u) - D 21 (u')D 12 (u)), 



where 



/l 



V 



a b 
b a 



x 



b{u, v) 



V 

b + (u, v)b~(u, v) 



/I 



V 



c d_ 



1/ 



, a(u, v ) = 1 — b(u, v) 



a (u,v)c + (u,v) 

By using these relations, we operate the transfer matrix (B.9) to the eigenstate 
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by (4.2). Here show the elements separately 

£>22(«)|$n(Ai,--- ,A„)> 

= (a+(u))in ^FM l$n(Al '''' ,An)> 

+ S^ =1 (a+(A,)) L |^ 1 2 1 KA J ;{A,})) 

+ E^EtJffxK A;, A,0(a + (A,)a + (A,0) L |*® i(«, A,-, A,; {A,})), (B.22) 

£>n(u)|$n(Ai,--- ,A„)) 
= (c + ( M )) L n^ 1 (-^^4)|$ n (A 1 , • • • , A n )) 

C yU, Aj) 

+ E^(6" (A,)) L A (1) (^ = K {A,})l*Si(«, A,; {A,})) 
+ E^Ci^K A,, A,)(6-(A0r (A,)) L 

x AW(« = A„ {A fe })A«(« = A,, {\ k })\^Uu, Xj, \r, {X k })), (B.23) 
^aA aa (m)|$„(Ai,-- - ,A n )) 

= (6- H) L n« =1 (-f^4)A«( M , {A i })|$ n (A 1) • • • , A n )) 

+ E? =1 (r (A,-))^* = A„ {X^Uu, A,; {A,})) 
+ ^ =1 (a + (X j )) L \^_ 1 (u,X f ,{X l })) 

+ E] =2 Ei-lH 3 (u, X t , XjMM, Xj) - a(X h A,-))(a + (A/)r (A,)) L 

x A( 1 )( M = A„{A fe })|< 3 2 1 ( M ,A„A,;{A fe })) 

+ £? =2 £Ci#4(«, A,, A,)(6(A„ AO - a(A„ A,))(a + (A j )r(A,)) L 

x M 1 \u = X l ,{X k })\¥^ 1 (u,X„X l ;{X k })}, (B.24) 
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where 



l*„-i(M*{M)> - -^n^.^^) 

® $ n "li(Ai, • • • , \j, ■ ■ ■ x n ).df\x f , {x k })3\o), 
a,; {A/})) = ^^n^.(-^4|)[f.(5H 1)] 

c + {u, Xj) ^ J b {X k ,Xj) 
<g> $„-i(Ai, • • • , Xj, ■ ■ ■ , X n ).6f\\ j ; {A fc }).#|0>, 

Aj ' A " {Xk})) ~ U ^> l b-(X k ,Xj) b~(\ k ,\i) 2l() * 

®4- 2 (Ai,-" ,A,v ,A^4f (A,,A j; {Aj).:?|0), (B.25) 

6f\X 3 ; {A,}) = na ^ (Afc '^ f fc , fc+1 (A fc , A,), (B.26) 

C^f (A;, Aj-; {X k }) = njjz\ _J J \r k +i,k+^k, Xj) 

x n£ = } +1 — — - k ,k+i(Xk, Aj)n^~\— y-^— ^-r fc,fc+i(A fe , A^), (B.27) 
a^(Ajfc, AjJ a^l^Afc, Aj) 



Hl( X yz)= a+ (y^ x ) e ( z ^ x )f(y^ x ) i d + (y,x)f(y,z) 
b-(y,x)b-(z,x)c + (y,x) c + (y,x)c + (y, z)' 
, . d + (x,y)f(y,z) e(x, y)b + (x, z) 

H2(x, y, z) = 1 '7 j + 7' ; - + yy + y y B.29 

ff3(x, y, z) = /(^y^^y) _ a + ^y)e{x,z)f{x y) 
c+(x,y)b {x,y)b {y,z) b {x,y)b (x,z)c + {x,y) 

rr A( \ f{x,y)e{x,y)e{y,z) e(y,x)e(x,y)f(y,z) 

n^\x,y,z) — z 



c+(x,y)b (x,y)b (y,z) b+(y,x)b (x,y)c+(y,z) 
a~(x,y)f(x,z)e(x,y) 



-(l + a(x,y)), (B.31) 



b (x,y)c+(x,z)b+(x,y) 
using the identities, 

4(Ai, • • • , A n ) = $ n (X j: Ai, • • • , Xj, ■ ■ ■ , X n ).6f\Xf, {X k }), (B.32) 

4(Ai, • • • , A n ) = ® n {X h Xj, Ai, • • • , A,, • • • , \, ■ ■ ■ , X^.dffiXi, A,; {X k }), (B.33) 
Hl(x,y,z) + H2(x,y,z) 

= (b(y, z) - a(y, z))H3(x, y, z) + (b(z, y) - a(z, y))HA(x, y, z). (B.34) 

The first term of r.h.s in (B.24) is obtained by using r-structure of the (B.ll) and the 
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following calculation. jF in (4.2) is fixed as follows, 

tT a tf\u,{\ l }).$ = A.W(u,{\ l })§, 

7^(u,{X l })=^(u,X n )---^(u,X 1 ), 

( 1 



5(w, Aj) a(w, Aj) 
d(u, \j) b{u, \j) 
\0 1/ 



(B.35) 



This L matrix XL*- 1 - 1 is a permutation of the r-matrix. After the reparametrization, 



a(A,/i) = 



U/2 



H = 



fx - X-U/2 
z-{n)-l/z-{n) , C/ 



, 6(A,/i) 



ji — \ 



p,- X-U/2' 



+ — , z (yu) = — exp(2/i(/i))/ tan^u, 



(B.36) 



we see that this is equivalent to the problem of diagonalizating the transfer matrix for the 
XXX model and we know the solution, 



aW(u, {a,}, M) = n^=-^— + u] =1 l(u, x^-J—. 



where I^ =1 b(^h \) = W&i 



b(/xi,/x k ) 



b(^ k ,/ii) 

Now the eigenvalue of the transfer matrix is obtained, 



(B.37) 
(B.38) 



AMM) = (a + («)) L II? =1 



^ a+ ^ u ) + ( c +(„))^n» =1 (-^M 

c (it, Xj) 



-(b-(u)) L UU( 



a~(u, Xj) 



)A«K{A,}). 



b (u,\j) 

And the unwanted terms in (B.22), (B.23) and (B.24) give the another constraint, 

,a+(A. 



(B.39) 



'6-(Ai) 



= A«( M = A„{M). 



(B.40) 
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